ST   ; TL99999;
V    1; W    0;
YA28671; YB    10; YC  999;
RESTART; J41; J19;

PROGRAM;

(P
 ABCSJL26.016
 Coset enumeration six cols fourth version.
 Original program by John Leech, 1965.
 This text has been formatted by WF to better display the program's structure.);

V0 = B1212121212121212; V1 = B0702/11;

30;(reader);
      SET 2; SET 5; OUT; = C15; SET AY0; = I15;
   (punch);
      SET 1; SET 5; OUT; = C14;
      C14; C14 TO Q13; SET 300; = M13; PGAP Q13;

19;(begin);
      SET AYB1;
   JS14; (read inverses);
   J 20;
   (all jobs done);
      PGAP Q13; C 13; SET 6; OUT; ZERO; OUT;

20;
   SET 1; DUP; = RM5; =C3; ZERO; DUP; =YA2; =YA3;

21;(read generators or relations);
      SET AYC0;
   JS14;
      ZERO; =RM4;

22;(next coset);
      M5; M4;
   J27=; (table complete);
      M+I4; ERASE; ZERO; =RM6;

23;(next relation);
      M4 TO Q2; YC0M6;
   J26=Z; (no more relations);

24;(next aR);
      YC0M6Q; DUP;
   J31=Z; (end of relation);
      =M1;
   JS1;
   J25; (no coset defined);
      =Q2; J24;

25;(define new coset);
      M5; SET 14335;
   J27=; (store full);
      ERASE; M+I5; M5;
   JS5; (insert forward reference);
      M5; DUP; +; =M3; ZERO; DUP; =YA0M3; =YA1M3;
      M2; M5 TO Q2; YB0M1; =M1;
   JS5; (insert backward reference);
   J24;

26;
   J22C3Z;
      DC3;
   J21; (subgroup defined, start of main work);

27;(output);
   JS99;
      (print number of cosets); ZERO; DUP; DUP; =RM2;

28;
      ERASE; ERASE; E0; SHL24;
   J19=Z; (test if table required);
      M2; M5; -;
   J19=Z; (table finished);
      M+I2; M2;
   JS99;
      M2; DUP; +; =M3; YA1M3; YA0M3;

29;
      DUPD; OR;
   J28=Z;
      ZERO; SHLD 16;
   JS100;
      SHL -16; SHLD 16;
   J29;

41;
      M5; J27;

1;(fetch);
      M2; DUP; +; =M3; M1;
   =LINK;
      YA1M3; YA0M3;

*2;
   EXIT AR2;
*  J3;
*  SHL 16; J3;
*  SHL 32; J3;
*  ERASE; J8;
*  ERASE; SHL 16; J8;
*  ERASE; SHL 32; J8;

3;
      REV; ERASE;
8;
      SHL -32; DUP;
   J4=Z;
   EXIT 2; (defined);

4;
      ERASE;
   EXIT 1; (undefined);

5;(store);
      M2; DUP; +; =M3; M1;
   =LINK;
      SHL 32;

*6;
   EXIT AR6;
*  J7;
*  SHL -16; J7;
*  SHL -32; J7;
*  J9;
*  SHL -16; J9;
*  SHL -32;

9;
      YA1M3; +; =YA1M3;
   EXIT 1;

7;
      YA0M3; +; =YA0M3;
   EXIT 1;

10;(delete);
      M2; DUP; +; =M3; M1;
   =LINK;
      ZERO; NOT; SHL 32;

*11;
   EXIT AR11;
*  J12;
*  SHL -16; J12;
*  SHL -32; J12;
*  J13;
*  SHL -16; J13;
*  SHL -32;

13;
      NOT; YA1M3; AND; =YA1M3;
   EXIT 1;

12;
      NOT; YA0M3; AND; =YA0M3;
   EXIT 1;

14;(input);
      =RM6; ZERO; =RM7; SET 125; I15; +; =M15; IM15 TO Q14;
      PREQ15; PWEQ14; Y0M7Q; DUP;

15;
      ERASE; ZERO; SHLD 6; SET 32; -; DUP;
   J15<Z; (ignore character);
      SET 29;
   J16 !=;
   EXIT 2; (all jobs done);

16;
      =M0M6Q; DUP;
   J17!=Z;
      ERASE; Y0M7Q;

17;
      ZERO; SHLD 6; SET 61;
   J18=; (end of sequence);
      SET 32; -; DUP;
   J16>Z;
      ERASE; ZERO; J16;

18;
      ERASE; ERASE; ZERO; DUP; =M0M6Q; =M0M6;
   EXIT 1;

99;(new line);
      ZERO; NOT; =TR;

100;(output integer);
      SET 1; =Q1; V0; REV; FRB; (number converted);

101;(add 16 to significant digit characters);
      Q1; -; DUP;
   J102<Z; (end of digits);
      Q1; +; Q1; SHA 6; DUP; =Q1; SHA -2; +;
   J101;

102;
      Q1; +;
   J103NTR;
      V1; +; (case normal, new line);

103;
      =Y0; I15; DUP; =I14; =M14; PWQ14;
   EXIT 1;

31;(look for coincidence);
     ERASE; M2; M4; -;
  J23=Z; (no concidence);
     M2; M4; MAX; SHL 24; +; =Y0;
     ZERO; DUP; =Y1; =Q7;
  JS50; (coincidence subroutine);
     M7; =RM8; ZERO; =RM12;

32;
      M8; DUP; +; =M3; YA0M3; YA1M3; OR;
   J35=Z;
      M8; SHL 24; M8; +; M12; -; =Y0;
   JS50; (use coincidence subroutine to close up);

36;
      M8; M4;
   J33!=;
      DUP; M12; -; =M4; (adjust current coset);

33;
      M5;
   J34!=;
      M12; -; =M5;
   J23; (adjust last coset defined);

34;  ERASE; M+I8; J32; (advance to next line);

35;
      M+I12;
   J36; (advance count of blank lines);

49;(vide infra);
      ERASE; M10 TO Q11;
   J53;

50;(coincidence subroutine);
      ZERO; Y0; SHLD -24; DUP; =Q9; Q7;
   J51=Z;
      Q7; MAX; ERASE;

51;
      =Q7; SHL -24; =Q10; (cosets in QQ9,10); ZERO; =RM1;

52;(next column);
      M+I1; M1; SET 7; -;
   J63=Z; (was last column);
      M9 TO Q2;
   JS1; (fetch bR);
   J52; (bR undefined);
      M9;
   J49=; (bR = b);
      =Q2; YB0M1; =M1;
   JS10; (delete bRR*);
      YB0M1; =M1; M2 TO Q11;

53;
      M10 TO Q2;
   JS1; (fetch aR);
   J59; (aR undefined);
      M9;
   J54!=;
      ERASE; M10; (if aR = b, replace by a);

54;(aR in N1);
      DUP; DUP; =M2; M11; ZERO; =RM11;

55;
      REV;
   J65=; (no new coincidences);
      MAX;

56;
      DUP; Y1M11Q; SHL -24; -;
   J56<Z; (search list);
      DUP; Y0M11; SHL -24; -;
   J57>Z; (new coincidence);
      ERASE; DUP; Y0M11; SHL 24; SHL -24;
   J55;

57;
      SHL 24; +;

58;
      Y0M11; REV; =Y0M11Q; DUP;
   J58!=Z; (move list down);
     =Y0M11; J60;

59;
      M11;
   JS5; (set aR = old bR);
      M11 TO Q2;

60;
      YB0M1; =M1;
   JS1; (fetch aRR*);
   J61;
   J62!=Z;

61;
      M10;
   JS5; (If aRR* undefined, set = a);

62;
      YB0M1; =M1; J52;

63;
      ZERO; DUP; M9; DUP; +; =M3; =YA0M3; =YA1M3; (delete coset b); ZERO; =RM11;

64;
      Y1M11; DUP; =Y0M11Q;
   J64!=Z; (move list up);
      Y0;
   J50!=Z;
   EXIT 1; (finished when list empty);

65;
      ERASE; ERASE;
   J60;

FINISH;

